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Abstract 

We introduce a new, fully massive, Monte Carlo program to compute all four- 
fermion processes in e + e~ collisions, including Higgs production. We outline 
our strategy for the matrix element evaluation, the phase space generation 
and the implementation of the leading higher order effects, and show, where 
available, comparisons with existing results. 
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1 Introduction 



Since the beginning of LEP2 four-fermion Physics underwent an extensive 
development. From an experimental point of view, new processes and effects, 
not included in the 1996 LEP2 Yellow Report Jl|, became relevant. All codes 
used so far have to be upgraded and extended with new features in order to 
match the improved experimental needs, especially in view of the final LEP2 
analysis. 

Although the needed improvements of the various event generators in use 
at LEP2 may differ from one code to another there often are three obvious 
steps to be taken, namely 

• Including fermion masses, instead of neglecting them. 

• Taking into account scale-dependent corrections and higher order con- 
tributions related to unstable particles. 

• Improving the treatment of the QED radiation. 

In the first place, fermion masses are relevant both for Higgs produc- 
tion and for single- W @ or 7 7 dominated processes || , when electrons are 
scattered in the very forward region. 

Secondly, the correct scales have to be taken into account for processes 
dominated by quasi-real t-channel photons (single- W production) or s-channel 
photons with low virtuality (Z7* processes), where tools adequate for study- 
ing high energy dominated processes, such as W + W~ production, fail in 
describing the data. 

Finally, it would be desirable to generate a realistic non-vanishing p t 
distribution for the radiated photons. 

These improvements should be made while, at the same time, the Monte 
Carlo program should be kept general enough to deal with all processes in 
all possible kinematical configurations. It is this requirement of having all 
processes and all kinematical configurations in one program, which is the 
challenge. 

This paper paves the way to replace an existing code, EXCALIBUR [|j], 
which describes all final states, by a new FORTRAN program, NEXTCALIBUR, 
which contains the above list of improvements. The demand of massive 
fermions leads us to a new type of matrix element calculation and severe 
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requirements on the phase space generation. This means in practice a new 
program, since the old program made essential use of massless fermions. 

The structure of the work is as follows. In section 2, we describe the 
strategy for the massive matrix element evaluation and for the Monte Carlo 
integration in NEXTCALIBUR. In section 3 we present the treatment of the 
leading higher order contributions, mainly QED radiation and running of 
olqed- Subsequently, in section 4, we show numbers and comparisons with 
other results in the literature. Finally, the last section is devoted to our 
conclusions. 

2 The strategy 

We solved the problem of the complexity in the massive matrix element 
evaluation by using HELAC, a FORTRAN package for helicity amplitudes com- 
putation based on Dyson- Schwinger equations, as described in ref. ||, to 
which we refer for more details. Here we only point out that, in the case of 
massless particles, the algorithm is simplified from the beginning, by avoiding 
the computation of helicity amplitudes that are known to be zero. A very 
fast computation of the cross section can therefore be obtained. In table 1 
we report comparisons in speed between NEXTCALIBUR (fully massive) and 
EXCALIBUR (massless) for three processes. 



Process 


NEXTCALIBUR 


EXCALIBUR 


e~ v e v e e + 


34.6 (16) 


12.3 (8) 


e~ v e u d 


35.0 (32) 


5.4 (3) 


fi + \l~ b b 


62.0 (64) 


6.1 (8) 



Table 1: CPU time comparison between NEXTCALIBUR and EXCALIBUR in 

seconds/1000 points. In parenthesis we show the number of non vanishing 
helicity configurations. The speed difference is mainly due to the increased 
number of contributing helicities. 

As for the phase space generation, we used the same multi-channel self- 
adjusting approach used in EXCALIBUR. Namely, we wrote a set of kinematical 
channels, each of them taking into account a different peaking structure of 
the integrand. The difference with the EXCALIBUR channels is that now all 
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Fzg 1: Higgs production (a) and double multi-peripheral (b) kinematical chan- 
nels in NEXTCALIBUR. 



fermion masses are taken into account. Furthermore the leading kinematical 
structures for Higgs production have been added (see fig. la). 

More details can be found elsewhere |J. We just mention here that 
the most complicated channel is the double multi-peripheral configuration 
given in fig. lb, because three t-channel like singularities have to be mapped 
simultaneously: two of them due to quasi-real photons, and the third one 
induced by the fermion propagator. 

3 Implementing higher order contributions 

A first source of numerically important higher order contributions comes from 
the widths of the unstable bosons, that have to be included without breaking 
gauge invariance. Our approach is to use complex bosonic masses everywhere, 
also in the definition of the weak mixing angle. This obeys all relevant Ward 
Identities and has been shown to be a very good approximation even in 
the forward region for t-channel dominated processes ||. 

A second source of large higher order effects is the QED radiation. A 
very common solution is using the Structure Function formalism, namely 
a convolution of the Born cross section together with QED Initial State 
Radiators ||. Such a strategy is implemented in most of the programs used 
for the analysis of the LEP2 data 0] and accurately reproduces the inclusive 
four-fermion cross sections, at least for s-channel dominated processes. 

Recent studies have shown, by inspection with the soft limit of exact 
calculations, that the Structure Function formalism can still be used for t- 
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channel dominated processes, provided the scale q 2 of the radiators is chosen 
to be of the order of the virtuality of the exchanged i-channel photons. Our 
approach to the QED corrections is therefore using Initial State Structure 
Functions with a proper choice of q 2 for each of the two incoming legs [J 

As to the scale q 2 , s should be taken for s-channel dominated processes, 
while, when a process is dominated by small t exchanges and —t is much 
smaller than s, the scale is related to t. This is e.g. the case in small 
angle Bhabha scattering and the proper scale is chosen as the one which 
reproduces roughly the exact first order QED correction, which is known 
for Bhabha scattering. A similar procedure now also exists for the multi- 
peripheral two photon process [12[], since an exact first order calculation is 



also available |TB[. In these t-channel dominated processes it is important 
to know whether a cross-section with angular cuts is wanted, since then the 
t-related scale will increase and the QED corrections as well. When no exact 
first order calculations are available the scale occurring in the first order soft 
corrections is also used as guideline to determine q 2 [fL2| [l4j| . 

In NEXTCALIBUR the choice of the scale is performed automatically by the 
program, event by event, according to the selected final state (see Table 2). 



Final State 


q 2 - 


q\ 


No e ± 


s 


s 


1 e~ 


|t_| 


s 


1 e+ 


s 


\u\ 


1 e~ and 1 e + 


\t-\ 


\t+\ 


2 e~ and 2 e + 


min( i_ ) 


min(|t+|) 



Table 2: The choice of the QED scale in NEXTCALIBUR. q\ are the scales 
of the incoming while t± represent the t-channel invariants obtained by 
combining initial and final state momenta. When two combinations are 
possible, as in the last entry of the table, that one with the minimum value 
of \t\ is chosen, event by event. 

However, as mentioned before, the choice of the appropriate q 2 is not 
the only required improvement to the treatment of the QED radiation. The 

Presently, Final State radiation is neglected, but it can be included with analogous 
techniques. 
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increasing precision of the collected LEP2 data also requires the knowledge 
of the p t spectrum of the emitted photons. We therefore use a particular 
form of p t dependent Structure Functions JI5| , derived, at the first leading 
logarithmic order, for small values of p t . 
In practice, we replace the quantity 

q 2 1 

m ( o) by J 

in the strictly collinear Structure Function for the i th incoming particle, by 
explicitly generating c\ and c 2 , the cosines (in the laboratory frame) of the 
emitted photons with respect to the incoming particles. Once are gener- 
ated, together with the energy fractions Xi j2 and the azimuthal angles 0i j2 , 
the momenta of two ISR photons are known. The four-fermion event is then 
generated in the c.m.s. of the incoming particles after QED radiation, and 
then boosted back to the laboratory frame. 



We also take into account non leading terms with the substitution [16 



Q 2 1 o m e 1 



The above choice ensures that the residue of the soft-photon pole gets 

2 

proportional to ln(-^) — 1, after integration over q. The inclusive QED 
result is therefore exactly recovered, after integrating over the pt spectrum 
and, at the same time, the pattern of the photon radiation is exact for small 
p t values. Notice that the p t spectrum is controlled by the same scale q 2 
used in the strictly collinear Structure Functions, namely an s-channel scale 
for s-channel processes and a t-channel scale for t-channel dominated final 
states. The radiation pattern is therefore different in the two situations, as 
one naively expects. The presence of such a scale can also be thought as 
an extra handle to tune our Monte Carlo predictions to the data, to get a 
satisfactory description of the radiation. 

Another problem, in presence of low t-channel scales, is that high energy 
renormalization schemes, such as the Gf scheme, fail in describing the data, 
because of the running of oiqed- 

A possible solution is the Fermion-Loop approach of refs. |l?|-||19fl, where 



all fermion corrections are consistently included by introducing running cou- 
plings g(s) and e(s) and re-summed bosonic propagators. 
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In presence of the WW'j vertex, the above ingredients are not sufficient 
to ensure gauge invariance, because loop mediated vertices have to be con- 
sistently included. On the contrary, when no WW'y vertex is present, the 
neutral gauge boson vertices, induced by the fermion loop contributions, are 



separately gauge invariant fl8 




Instead of explicitly including the loop vertices, we follow a Modified 
Fermion-Loop approach. Namely, we neglect the separately gauge invariant 
neutral boson vertices, and include only the part of the WW'j loop function 
necessary to renormalize the bare WW'y vertex and to insure the U(l) gauge 
invariance []. Our procedure is as follows. Besides running couplings, we use 
bosonic propagators 

Pr(s)=(s-M 2 z (s)Y 1 (g» u - 
with running boson masses defined as 

M l{s) = fi w 9 2( ^ S \ - g 2 (s)[T w (s) - T w (fx w )] 

Tw,z(s) are contributions due to the top quark, fj, W)Z the complex poles of the 
propagators (one can take, for instance, fi w>z = z — iT W)Z M WjZ ) and 

400 = 44, <$(*) = 1 - sj(s) . 
9 \ s ) 

The explicit form of the running functions e 2 (s), g 2 {s) and Tw,z{s) can be 
found in ref. |T8|. 



The leading contributions are in the real part of the running couplings 
therefore we take only the real part of them. This also means that one can 

2 When using NEXTCALIBUR in this running coupling mode, U(l) is preserved but St/ (2) 
is, in general, violated. The numerical effects of this are expected to be small at LEP2 
energies. The implementation in the code of the fully SU(2) x t/(l) gauge invariant solution 
of ref. pfJ is under way. 
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replace, in the above formulae, g 2 ([i WtZ ) — > g 2 (M 2 z ), c^(// z ) — > Cg(M 2 ) and 
also TV^/i^) -> Tw^M^J. 

When the WW'j coupling is present, we introduce, in addition, the fol- 
lowing effective three gauge boson vertex 




with s = p 2 , s" 



P 



+ ■ 



p 2 _ and 



V, 



Hup 



+ 



'V 



9»v(P ~ P+) P + 9v P {p+ ~ P-)p (l + S v ) + g P p{p- - p) v 



(p + -p_). 



s — s~ 



1 P-uP-c 



g(s + )g(s )(s - s , 

+ - g(s-)] [s-g(s + ] + 



~g 2 (s+)g 2 (s-){T w (s~)-T w (s + )] 



(1) 



It is the easy to see that, with the above choice for V^ up , the U(l) gauge 
invariance - namely current conservation - is preserved, even in presence of 
complex masses and running couplings, also with massive final state fermions. 

From eq. (1), one deduces at least two effective ways to preserve U(l). 
One can either compute g(s) at a fixed scale (for example always with s = 
Myy), while keeping only the running of e(s), or let all the couplings run at 
the proper scale. With the first choice the modification of the three gauge 
boson vertex is kept minimal (but the leading running effects included). With 
the second choice everything runs, but a heavier modification of the Feynman 
rules is required. At this point one should not forget that our approach is an 
effective one, the goodness of which can be judged only by comparing with the 
exact calculation of ref. [fTTp . We found that the second choice gives a better 
agreement for leptonic single- W final states, while the first one is closer to the 
exact result in the hadronic case, which is phenomenologically more relevant. 
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Therefore, we adopted this first option as our default implementation in 
NEXTCALIBUR. 

We want to stress once more that the outlined solution is flexible enough 
to deal with any four-fermion final state, whenever small scales dominate. 
For example, once the given formulae are implemented in the Monte Carlo, 
the correct running of ciqed is taken into account also for s-channel processes 
as Z7* production. 

4 Numerical results 

In tables 3 and 4, we show the total cross sections for the processes e + e" — > 
e + e~fi + fi~ and e + e~ — > e + e~e + e~ . Where available, we compare our predic- 
tions with the QED numbers published in ref. || . 





BDK 


NEXTCALIBUR 


20 


98.9 ± 0.6 


99.20 ± 0.98 


35 


131.4 ± 2.2 


131.03 ± 0.88 


50 


154.4 ± 0.9 


152.33 ± 0.83 


100 


205.9 ± 1.2 


204.17 ± 1.73 


200 




263.50 ± 1.31 


200 (all) 




265.58 ± 1.44 



Table 3: <Jtot (in nb) for the process e + e — > e + e [i + fi . Only QED diagrams, 
except in the last entry. 

NEXTCALIBUR contains all electroweak diagrams, and can therefore be used to 
compute the electroweak background to the above 7 7 processes. By looking 
at the last entry of the tables, the latter is found to be less than 1 % at LEP2 
energies, at least for totally inclusive quantities. 

All numbers have been produced at the Born level, but ISR and running 
ctQED can be included as described in the previous section. Here we mainly 
want to demonstrate the ability of the program to cover all phase-space 
regions, without loosing efficiency in the event generation. 

At the moment, in order to get the necessary numerical accuracy, we run 
the program in quadruple precision. NEXTCALIBUR has been written in such a 
way that switching from double to quadruple precision simply implies adding 
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BDK 


NEXTCALIBUR 


20 


0.920 ± .011 


0.905 ± .011 


35 


1.070 ± .015 


1.079 ± .014 


50 


1.233 ± .018 


1.214 ± .016 


100 


1.459 ± .025 


1.485 ± .020 


200 




1.776 ± .019 


200 (all) 




1.787 ± .030 



Table J^: a tot (in nb xl0 7 J for the process e + e — > e + e e + e . Only QED 
diagrams, except in the last entry. 



a flag at the compilation time. However, this option is really necessary only 
when two or more electrons are allowed in the very forward direction. For 
all the other kinematical configurations, with at most one electron lost in 
the beam pipe, double precision is sufficient. A version of the program using 
double precision in all possible situations is currently under study. 

In tables 5, 6 and 7 we show, as a second example, single-jy numbers 
produced with our Modified Fermion-Loop approach, as discussed in the 
previous section. Comparisons are made with the exact Fermion-Loop calcu- 
lation of ref . |l7j . The results of the complete Fermion-Loop are reproduced 
at 2% accuracy for both leptonic and hadronic single-H^ final states. 

It should also be noted that, when neglecting Fermion-Loop corrections, 
one can directly compare NEXTCALIBUR with other massive Monte Carlo's 
and one finds excellent agreement for single- W production in the whole phase 
space 
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In Figs. 2 and 3, we show, as an illustrative example, two distributions 
for the most energetic radiated photon in the processes e + e~ — > fi^v^ud^), 
as predicted by NEXTCALIBUR. Only ISR photons are taken into account, 
y/s = 200 GeV, |cos^| < 0.985, E^ > 5 GeV and M(ud) > 10 GeV. 
With the same set of events, by using as separation cuts for the emitted 
photons E 7 > 1 GeV and |cos# 7 | < 0.985, we found the following values 
for the total, the non-radiative, the single-radiative (1 generated 7) and the 
double-radiative (2 generated 7's) cross sections: 

a tot = 0.61727 ± 0.0059 pb a n _ rad = 0.57819 ± 0.0058 pb 
o-s-rad = 0.03854 ± 0.0016 pb a d _ rad = 0.00054 ± 0.0002 pb 
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da /d6 e 


MFL 


EFL 


MFL/EFL - 1 (percent) 


0.0° -=-0.1° 


0.45062(70) 


0.44784 


+0.62 


0.1° -=-0.2° 


0.06636(28) 


0.06605 


+0.47 


0.2° -=-0.3° 


0.03848(21) 


0.03860 


-0.31 


0.3° -=-0.4° 


0.02726(18) 


0.02736 


-0.37 


d tot 


83.26(9) 


83.28(6) 


-0.02 



Table 5: da/d6 e [pb / 'degrees] and a to t [ft] for the process e + e — > e v e ud. 
The first column is our Modified Fermion-Loop, the second one is the exact 
Fermion-Loop of ref. |T/|/. i/i = 183 GeV, |cos# e | > 0.997, M(ud) > 45 
GeV. QED radiation not included. The number in parenthesis shows, when 
available, the integration error on the last digits. 



da / d6 e 


MFL 


EFL 


MFL/EFL - 1 (percent) 


0.0° ^0.1° 


0.13218(26) 


0.13448 


-1.7 


0.1° ^0.2° 


0.01997(10) 


0.02031 


-1.7 


0.2° ^0.3° 


0.01171(8) 


0.01194 


-1.9 


0.3° -=-0.4° 


0.00838(6) 


0.00851 


-1.5 


CTtat 


25.01(3) 


25.53(4) 


-2.0 



Table 6: da /d6 e [pb/degrees] and a to t [fb] for the process e + e — > e z/gZ/^/i" 1 ". 
The first column is our Modified Fermion-Loop, the second one is the exact 
Fermion-Loop of ref. y/s = 183 GeV, | cos 6> e | > 0.997, |cos^| < 0.95 

and > 15 GeV. QED radiation not included. The number in parenthesis 
shows, when available, the integration error on the last digits. 
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da j d8 e 


MFL 


EFL 


MFL/EFL - 1 (percent) 


0.0° -=-0.1° 


0.62694(98) 


0.62357 


+0.54 


0.1° ^0.2° 


0.08850(38) 


0.08798 


+0.59 


0.2° -=-0.3° 


0.05100(30) 


0.05141 


-0.80 


0.3° -=-0.4° 


0.03672(25) 


0.03646 


+0.71 


O~tot 


113.73 (13) 


113.67(8) 


+0.05 



Table 1: da/d6 e [pb / 'degrees] and a to t [fb] for the process e + e — > e v e ud. 
The first column is our Modified Fermion-Loop, the second one is the exact 
Fermion-Loop of ref. ^j. yfe = 200 GeV, \ cos6 e \ > 0.997, M(ud) > 45 
GeV. QED radiation not included. The number in parenthesis shows, when 
available, the integration error on the last digits. 



E — i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 3 



1 da 



a d cos 9~, 



10 1 



10 c 



io- 1 - 




I I I I I I I I I I I I I I I I I I 1_ 

1 1 

cos6L 



Fig 2: cos# 7 distribution (with respect to the incoming e + ) for the most 
energetic photon in the process e + e~ — > fi~ u^ud^) . 
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1 da 

a dE~< 



10° 



10- 1 : 



10"* - 



10" 3 - 



io~ 4 - 




10 



20 30 
£ 7 [GeV] 



40 



Fig 3: E 1 distribution for the most energetic photon in the process 



e T e" 



Finally, in tables 8 and 9, we show comparisons with the Higgs cross 
sections published in ref. |21|| , by choosing WPHACT as a benchmark 
program. We devoted special care to implement exactly the same input 



parameters of ref. ||21|| . For completeness we list them here: 

• Standard LEP2 input parameter set (see ref. H). 

• Massless fermions everywhere, except in the Higgs coupling to the b. 

• Running widths in the bosonic propagators. 
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rein ^vjev ) 


65 


90 


115 


Final state 


bb 


WPHACT 


32.7141(68) 


1.59946(64) 


1.05953(56) 


NEXTCALIBUR 


32.691(19) 


1.5999(18) 


1.0588(11) 


Final state 




WPHACT 


64.238(14) 


2.3661(10) 


1.29237(82) 


NEXTCALIBUR 


64.256(35) 


2.3651(27) 


1.2910(14) 


Final state 


v e v e bb 


WPHACT 


71.694(27) 


5.0996(23) 


1.08027(89) 


NEXTCALIBUR 


71.778(84) 


5.086(12) 


1.0776(13) 



Table 8: Cross sections, in fb, for Higgs production at yfs= 175 GeV. 
See text for the input parameters. 

• m T = 1.777 GeV, m b = 2.9 GeV and m c = 0.75 GeV. 

• No ISR, no QCD corrections but all background diagrams included. 

• M z - 25 GeV < m a < M z + 25 GeV and m bl > 50 GeV. 

It is worth mentioning explicitly that NEXTCALIBUR can consistently include 
fermion masses everywhere, and that they have been neglected in the pre- 
sented numbers just for the sake of comparison. 

5 Conclusions 

We introduced NEXTCALIBUR, a new Monte Carlo program to study four- 
fermion processes in e + e~ collisions. We outlined our strategy for includ- 
ing Higgs, fermion masses and leading higher order effects, without loosing 
efficiency in the event generation. The program is meant to upgrade the 
performances of an already existing code 0]. 

We concentrated mainly on QED and scale-dependent corrections, with- 
out making any attempt to include genuine weak contributions. While the 
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rein ^vjev ) 


65 


90 


115 


Final state 


H + H bb 


WPHACT 


37.3990(64) 


24.4727(40) 


10.7027(24) 


NEXTCALIBUR 


37.394(21) 


24.471(14) 


10.7006(77) 


Final state 




WPHACT 


72.927(16) 


47.222(12) 


19.841(11) 


NEXTCALIBUR 


72.929(46) 


47.231(33) 


19.842(19) 


Final state 


v e v e bb 


WPHACT 


80.611(34) 


53.335(19) 


20.893(12) 


NEXTCALIBUR 


80.507(96) 


53.280(67) 


20.897(24) 



Table 9: Cross sections, in fb, for Higgs production at yfs= 192 GeV. 
See text for the input parameters. 



latter are certainly relevant for LEP2 precision measurements, such as aww 
and Mw, they do not seem to be necessary for all the other observables. In 
that respect NEXTCALIBUR represents a solid tool for the final analysis of the 
LEP2 four-fermion data and for studying e + e~ Physics in general. 

In the near future, two big improvements of the program are foreseen. 
First of all the inclusion in the matrix element of anomalous couplings. Sec- 
ondly, the implementation of the formalism of ref. [20 to incorporate run- 
ning couplings and finite boson widths effects without breaking SU(2) x {7(1) 
gauge invariance. 

Notice that, given our computational strategy, the only needed modifi- 
cation is the insertion of additional Feynamn rules in the matrix element, 
all the rest remaining the same. Therefore, we do not expect difficulties of 
principle in the actual implementation. 
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